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Abstract. Four members of the new class of pulsating sdB stars are 
analysed from Keck HIRES spectra using NLTE and LTE model at- 
mospheres. Atmospheric parameters (T e ff, log g, log(He/H)), metal 
abundances and rotational velocities are determined. Balmer line fit- 
ting is found to be consistent with the helium ionization equilibrium 
for PG 1605+072 but not so for PG 1219+534 indicating that sys- 
tematic errors in the model atmosphere analysis of the latter have 
been underestimated previously. All stars are found to be helium de- 
ficient probably due to diffusion. The metals are also depleted with 
the notable exception of iron which is solar to within error limits 
in all four stars, confirming predictions from diffusion calculations of 
Charpinet et al. (1997). While three of them are slow rotator's (vsini 
< lOkm/s), PG 1605+072 displays considerable rotation (vsini = 
39km/s, P<8.7h) and is predicted to evolve into an unusually fast 
rotating white dwarf. This nicely confirms a prediction by Kawaler 
(1999) who deduced a rotation velocity of 130km/s from the power 
spectrum of the pulsations which implies a low inclination angle of 
the rotation axis. 
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1. INTRODUCTION 

It is now well established that the hot subluminous B stars can 
be identified with models of the extreme Horizontal Branch (EHB) 
stars (Heber, 1986, Saffer et al., 1994). 

Recently, several sdB stars have been found to be pulsating 
(termed EC14026 stars after the prototype, see O'Donoghue et al. 
1999 for a review), defining a new instability strip in the HR-diagram. 
The study of these pulsator's offers the possibility of applying the 
tools of asteroseismology to investigate the structure of sdB stars. 
The existence of pulsating sdB stars was predicted by Charpinet et 
al. (1996), who uncovered an efficient driving mechanism due to an 
opacity bump associated with iron ionization in EHB models. How- 
ever, in order to drive the pulsations, iron needed to be enhanced 
in the appropriate subphotospheric layers, possibly due to diffusion. 
Subsequently, Charpinet et al. (1997) confirmed this assumption 
by detailed diffusion calculations. Even more encouraging was the 
agreement of the observed and predicted instability strip. 

Thirteen pulsating sdB stars are well-studied photometrically 
(O'Donoghue et al. 1999). A precise knowledge of effective temper- 
ature, gravity, element abundances and rotation is a prerequisite for 
the asteroseismological investigation. 

We selected four EC14026 stars for a detailed quantitative spec- 
tral analysis: PG 1605+072 was chosen because it has the lowest 
gravity and, therefore, has probably already evolved beyond the ex- 
treme horizontal branch phase. It also displays the richest frequency 
spectrum amongst the EC 14026 stars (>50 periods have been iden- 
tified, Kilkenny al. 1999). Recently, Kawaler (1999) predicted from 
his modelling of the pulsations that PG 1605+072 should be rotat- 
ing. PG 1219+534 was chosen because it has the shortest pulsation 
periods and has a helium abundance larger than most other sdB stars 
(O'Donoghue et al., 1999). For Feige48 and KPD 2109+4401 only 4 
rsp. 5 frequencies have been found so far. Feige48 is also the coolest 
of all EC 14026 stars known. 

2. OBSERVATIONS 

High resolution optical spectra of the four pulsating sdB stars 
were obtained with the HIRES echelle spectrograph (Vogt et al. 
1994) on the Keck I telescope on July 20, 1998 using the blue cross 
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disperser to cover the full wavelength region between 3700A and 
5200A at a resolution 0.09A. 

The spectra are integrated over one pulsation cycle or more since 
the exposure times (600-900s) were long compared to the pulsational 
periods. 

The standard data reduction as described by Zuckerman & Reid 
(1998) resulted in spectral orders that have a somewhat wavy con- 
tinuum. In order to remove the waviness we used the spectrum of 
H1504+65 (a very hot pre-white dwarf devoid of hydrogen and he- 
lium) which was observed in the same night. Its spectrum has only 
few weak lines of highly ionized metals in the blue (3600-4480A) 
where the strong Balmer lines are found in the sdB stars. Therefore 
we normalized individual spectral orders 1 to 20 (3600-4480A) of the 
sdB stars by dividing through the smoothed spectrum of H1504+65. 
The remaining orders were normalized by fitting the continuum with 
spline functions (interpolated for orders 26 and 27 which contain 
H/3). Judged from the match of line profiles in the overlapping parts 
of neighboring orders this procedure worked extremely well. At- 
mospheric parameters determined from individual Balmer lines are 
found to be consistent with each other except for H/3. Therefore, 
we excluded H/3 from the fit procedure. Moreover, the resulting T e fj 
and log g are also in excellent agreement with those from the fit of a 
low resolution spectrum of PG 1605+072 obtained at the ESO NTT. 
Details on the analysis of PG 1605+072 can be found in Heber et al. 
(1999a). 

3. ATMOSPHERIC PARAMETERS 

The simultaneous fitting of Balmer and He line profiles by a 
grid of synthetic spectra (see Saffer et al. 1994) has become the 
standard technique to determine the atmospheric parameters of sdB 
stars. The Balmer lines (H7 to H12), He I (4471 A, 4026A, 4922 A, 
4713A, 5016A, 5048A) and He II 4686A lines are fitted to derive all 
three parameters simultaneously. 

The analysis is based on grids of metal line blanketed LTE model 
atmospheres for solar metalicity and Kurucz' ATLAS6 Opacity Dis- 
tribution Functions (see Heber et al. 1999b). Synthetic spectra 
are calculated with Lemke's LINFOR program (see Moehler et al. 
1998). In addition, a grid of H-He line blanketed, metal free NLTE 
model atmospheres (Napiwotzki 1997), calculated with the ALI code 
of Werner & Dreizler (1999). 
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PG 1605+ 072 

The results (T eff =31 900K, log g=5.29, log(He/H)=-2.54) are in 
agreement with those from low resolution spectra analysed with sim- 
ilar models (Koen et al. 1998) as well as from our own low resolution 
spectrum for PG 1605+072. 

Four species are represented by two stages of ionization (He I 
and He II, C II and C III, N II and N III, Si III and Si IV). Since 
these line ratios are very temperature sensitive at the temperatures in 
question, we alternatively can derive T e fj and abundances by match- 
ing these ionization equilibria. Gravity is derived subsequently from 
the Balmer lines by keeping T e ff and log(He/H) fixed. These two 
steps are iterated until consistency is reached. C II is represented by 
the 4267A line only, which is known to give notoriously too low car- 
bon abundances. Indeed the carbon ionization equilibrium can not 
be matched at any reasonable T e ff. The ionization equilibria of He, 
N and Si require T e ff to be higher than from the Saffer procedure, 
i.e. 33 200K (He), 33 900K (N) and 32 800K (Si). 

Since this difference could be caused by NLTE effects, we re- 
peated the procedure for T e ff and log(He/H) using NLTE models. 
Alternatively, applying Saffer's procedure with the NLTE model grid 
(see Fig. 1) yields T e ff almost identical to that obtained with the 
LTE grid. Evaluating the He ionization equilibrium in NLTE, in- 
deed, results in T e ff being consistent with that from Saffer's proce- 
dure. We therefore conclude that the higher T e ff derived above from 
the ionization equilibrium in LTE is due to NLTE effects. 

However, a systematic difference in log g persists, the LTE values 
being higher by 0.06 - 0.08 dex than the NLTE results. Since its ori- 
gin is obscure, we finally adopted the averaged atmospheric param- 
eters: T eff =32 300±300K, log g=5.25±0.05, log(He/H)=-2.53±0.1. 
Helium is deficient by a factor of 30 as is typical for sdB stars. 

Feige48 

Since no He II line can be detected, the helium ionization equi- 
librium can not be evaluated. The procedure of Saffer et al. (1994) 
results in T eff =29400K, log g=5.51, log(He/H)=-2.90. Feige48 has 
the lowest helium abundance among our programme stars. 
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Fig. 1. Balmer and He line profile fits for PG 1605+072 of the HIRES 
spectrum from NLTE model atmospheres. 



KPD 2109+4401 

The helium ionization equilibrium and the Saffer et al. procedure 
give (averaged) parameters T ef f=31800K, log g=5.79, log(He/H) = 
-2.22 for KPD 2109+4401. 
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PG 1219+534 

Unlike for PG 1605+072, the helium ionization equilibrium and 
the Saffer et al. procedure give discrepant results: T e ff=33 200K, 
log g=5.93, log(He/H)=-1.60 (Saffer et al. procedure, Fig. 2) and 
T eff =35 200K, log g=6.03, log(He/H)=-1.41 (He ionization equilib- 
rium, Fig. 3). At the lower T e ff the Balmer lines are well matched 
through out the entire profile, whereas for He II 4686A there is a 
significant mismatch (see Fig. 2). At the higher T e ff He II 4686A is 
well reproduced, but the Balmer line cores are not reproduced at all 
(see Fig. 3). Despite of its high gravity PG 1219+534 has an unusu- 
ally high helium abundance, i.e. it is deficient by a factor of 2 to 5 
only. The line cores of He I 4026A and 447lA cannot be reproduced 
by either model. We conclude that our models do not describe the 
outermost layers of the atmosphere correctly where the cores of the 
Balmer and He I lines are formed. We point out that PG 1219+534 
has the highest helium abundance and the shortest pulsation periods, 
which might affect the outermost layers. 

4. ABUNDANCES 

Weak metal lines are present in the spectra of all programme 
stars. However, the number of detectable lines differs considerably. 
The largest number of metal lines is present in Feige48 (C, N, O, 
Ne, Mg, Si, Al, S and Fe) and PG 1605+072 (which lacks Al and S). 
In PG 1219+534 only N, S and Fe are detectable. 

The metal lines are sufficiently isolated to derive abundances 
from their equivalent widths except for the crowded region from 
4635A to 4660A in PG 1605+072 which we analyse by detailed spec- 
trum synthesis. Results are plotted in Figure 4. Upper limits are 
shown when no line of a species was detectable. Although several 
O lines are available in the spectra of PG 1605+072 and Feige48, it 
was impossible to determine the microturbulent velocity in the usual 
way, i.e. by minimizing the slope in a plot of the O abundances ver- 
sus equivalent widths, due to the lack of sufficiently strong lines. 
We adopted 5±5km/s which translates into small systematic abun- 
dance uncertainties of ±0.05dex for most ions. The analysis is done 
in LTE. A temperature uncertainty of AT e ff=1000K translates into 
abundance uncertainties of less than 0.1 dex. Hence systematic errors 
are smaller for most ions than the statistical errors. 
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Fig. 2. Balmer and He line profile fits for PG 1219+534 of the HIRES 
spectrum. Note the mismatch of the He II 4686A line profile and the cores 
of He I 4026A and 447lA. 

Like helium the metals are deficient with the notable exception 
of iron, which is solar to within the error limits. The high gravity 
stars (KPD 2109+4401 and PG 1219+534) have considerably lower 
O and Si abundances than the stars of somewhat lower gravity which 
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Fig. 3. He line profile fits for PG 1219+534 of the HIRES spectrum to 
determine T e fj and log(He/H) simultaneously, log g is adjusted to match 
the Balmer line wings. Note the mismatch of the cores of the Balmer lines 
and of He I 4026A and 447lA. 



point to the (selective) action of diffusion. It is, however, puzzling 
that iron is solar irrespective of the stellar gravity. UV spectroscopy 
is required to determine more precise iron abundances. We point 
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Fig. 4. Abundances of PG 1605+072 relative to the sun. Upper limits 
are denoted by arrows. 

out that a solar surface abundance is in perfect agreement with the 
diffusion calculations of Charpinet et al. (1997). 
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5. ROTATION VELOCITIES 

The spectral lines of PG 1605+072 are considerably broadened, 
which we attribute to stellar rotation and derive vsini = 39km/s, by 
fitting the strongest metal lines. In Figure 5 we compare a section 
of the spectrum of PG 1605+072 to that of the Feige48, which (like 
PG 1219+534 and KPD 2109+4401) are very sharp-lined (vsini<8- 
lOkm/s). 

Assuming a mass of 0.5 M the radius of R=0.28R© for 
PG 1605+072 follows from the gravity. Since sini cannot be con- 
strained the corresponding rotation period of PG 1605+072 must be 
smaller than 8.7h. PG 1605+072 displays the most complex power 
spectrum with more than 50 frequencies identifiable (Kilkenny et al. 
1999), 39 being bona fide normal pulsation frequencies. 

Usually rotation becomes manifest in the power spectrum by 
the characteristic splitting into equidistantly spaced multiplet com- 
ponents as is observed e.g. for the pre-white dwarf PG 1159—035 (ro- 
tation period: 1.4 d, Winget et al. 1991). Such multiplet's, however, 
have not been identified for PG 1605+072. Fast rotation introduces 
higher order terms that result in unequally spaced multiplet compo- 
nents. Recently, Kawaler (1999), was able to identify the five main 
peaks by considering mode trapping and rotational splitting. He pre- 
dicted that PG 1605+072 should be rapidly rotating (130 km/s). The 
measured v sin i= 39 km/s, hence, is a nice confirmation of Kawaler's 
prediction. Taken at face value a low inclination angle of 17 degrees 
results. 

Rotation is interesting also from the point of view of stellar evo- 
lution. PG 1605+072 is probably already in a post-EHB phase of 
evolution (Kilkenny et al. 1999) and will evolve directly into a white 
dwarf, i.e. will shrink from its present radius of 0.28 R© to about 
0.01 R©. Hence PG 1605+072 will end its life as an unusually fast 
rotating white dwarf if no loss of angular momentum occurs. Isolated 
white dwarfs, however, are known to be mostly very slow rotators 
(e.g. Heber et al. 1997, Koester et al. 1998). 

ACKNOWLEDGEMENT. 

U.H. gratefully acknowledges financial support by NATO ARW 
funds. 



Spectral analysis of pulsating sdB stars 



11 




4640 



Fig. 5. Fit of a section of the metal line spectrum of PG 1605+072, 
(bottom, v sini=39 km/s) compared to that of Feige48 (top, no rotation). 
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